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ABSTRACT 


In  the  synthesis  of  sampled-data  control  systems, 
the  criteria  of  performance  are  expressed  in  the  form 
of  time  domain  specifications,  such  as  overshoot,  peak 
time  and  damping  ratio.  Because  of  the  flexibility  in 
the  application  and  realizability  of  a  digital  controller, 
a  sampled-data  system  can  be  designed  to  exhibit  deadbeat 
response  whenever  a  specific  test  signal  such  as  a  step, 
ramp  or  parabolic  input  is  applied. 

If  however,  the  system  contains  nonlinear  elements, 
the  design  of  the  digital  controller  depends  not  only 
upon  the  type  of  test  signal  applied  but  upon  the  magnitude 
as  well. 

In  this  thesis  a  'feedforward’  method  is  used  to 
compensate  for  the  nonlinearity,  so  that  the  system 
exhibits  deadbeat  response  regardless  of  the  magnitude 
of  the  step  input. 

Most  references  on  the  subject  consider  systems  in 
which  the  nonlinear  element  can  be  isolated  from  the 
remainder  of  the  plant.  In  most  practical  situations 
this  is  not  the  case  and  the  nonlinearity  is  in  fact  an 
integral  part  of  the  plant.  This  thesis  considers  both 
cases  and  develops  two  methods  for  determining  the  type  of 
feedforward  compensation  required  in  each  case. 
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Although  perfect  deadbeat  response  was  not 
obtained  for  large  signals,  a  marked  improvement 
was  obtained  over  the  case  in  which  no  compensation 
was  used. 
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CHAPTER  ONE 


INTRODUCTION 

In  the  study  and  design  of  feedback  control 
systems,  performance  specifications  cannot  usually 
be  met  merely  by  adjusting  certain  parameters  of 
the  existing  system.  The  study  of  compensation 
techniques  then  becomes  of  great  importance. 

The  problem  which  is  considered  in  this  thesis 
is  that  of  compensating  a  nonlinear  sampled-data 
system  so  that  it  exhibits  deadbeat  response  for  a 
step  input.  A  linear  system  can  be  compensated  so 
that  the  output  exhibits  minimal  (or  deadbeat) 
response  regardless  of  the  magnitude  of  the  step 
input.  For  a  nonlinear  system,  however,  the  comp¬ 
ensation  required  is  dependent  upon  the  magnitude 
of  the  step  input. 

In  sampled-data  systems,  the  principle  of 
compensation  is  similar  to  that  for  continuous -data 
systems.  The  output  response  in  both  cases  is 
characterized  by  the  usual  terms  of  specifications 
such  as  per-cent  overshoot,  rise  time,  steady 
state  error  and  so  forth.  Because  of  the  sampling 
operation,  the  design  of  compensating  devices  for 
sampled-data  systems  is  usually  more  complicated 
than  for  continuous-data  systems.  However,  because 


V 


3H0  A3T3AH3 

.  '•  HOiTOuaojmii 

■ 

YlXjBiiSxj  :)cnnbo  anp.U  ^ofmrato;>  ec  t;  v  .v-c 

. 

. 

alaaxfr  .9±th  rrx  bo*  ,h-iaxioo  aj  rT?.;dw  ns  idoiq  &  r< 

..  -».d'i0O/f  /-,  %xi ::■&?■, Ci” HP*  ■  -  '•  •  r  : 

b  *;o^  o8f30<  lB9drjEcfa  adJ  '•  9  -  ipc ;  r *  -»^or<;a 

03  ba^j  *jor>  3d  nso 

(j  jsdbfidfe  10)  /  m  r  !  g  :  xo  ijjq^uo  Z:l 

CT93  3  !Xfo  lo  obi/*!  Igl  lo  -  iiblfcgV:;  S^OvBO* 

■ 

, 

. 

5o  .3/  :,  orJiq  r  ,  vne.  aye  ?-  —l  >c<  *  .  C.  .  ^ 

.^1  Irvjcufl  105  :i  a  3  Bill  OJ  /?;j.  •  3q: (GO  . 

ai  89BB3  ffrod  fri  3&ttoq?9*  b  ■  odl  .  am^ava 

! 

,bB.iJa  «.o(Tiii  9^'  ,oood.  iiVv*-  . .  >0 -toe  as  nous 
nnxlcrm  >xk)  o  ■.  ii'f.  ■  I  .d  n  -  o  .  b.  ,  ten  ,.  £-J  9 

10^  g33XV3b  grriosax* 
b90BO 

oa  >:  ■.  .  3V  9/ed  >3 .no  o.,iJ'~>*.  ■  :>  toi, 


1 


■a 


of  sampling,  the  design  procedure  is  usually  much 
more  versatile. 

The  method  of  compensation  which  will  be 
considered  in  this  thesis  is  accomplished  by 
introducing  digital  controllers  into  the  system. 

A  digital  controller  is  simply  a  pulsed-data  net¬ 
work,  of  which  the  transfer  function  can  be  expressed 

hy:  -i  -o 

D(z)  =  fc_+  alz_  +  a2z  +  •  •  •  +  amz _ 

b0  +  btz  +  b2z  +  ...  +  bnz'n 

where  n  is  larger  than  or  equal  to  m.  The  terms 

in  the  numerator  represent  the  inputs  at  the 

various  sampling  instants.  Similarily,  the  terms 

in  the  denominator  represent  the  values  of  the 

output  of  the  digital  controller  at  the  sampling 

instants.  In  general  n  is  equal  to  m  and  since 

the  input  and  output  values  are  valid  only  at 

the  specific  sampling  instants,  the  transfer 

function  can  be  written  as : 

D(s)  -  K0z°  +  K^z  ^  +  . ..  +  Y^z  n  . 

The  digital  controller  can  now  simply  be  regarded 
as  a  variable  gain  element,  in  which  the  gain  is 
constant  throughout  any  sampling  period. 

In  this  thesis,  the  following  general 
configuration  shall  be  considered.  (Fig.  1-1) 
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Fig *  1-1 

The  desired  output  will  be  deadbeat  response  and 
the  only  type  of  input  signal  which  will  be  con¬ 
sidered  is  a  step  input, 

1=1  THE  PLANT 

The  plant  is  a  simple  Type  1  system  with  a 
transfer  function,  G(s)  =  l/s(s+l).  In  addition, 
however,  it  will  contain  a  saturation  type  of 
nonlinearity.  Two  situations  will  be  considered0 
The  first  case,  in  which  the  nonlinearity  is  iso¬ 
lated  from  the  plant,  is  the  case  which  most  authors 
consider  in  examples  of  nonlinear  compensation.  The 
approach  is  usually  quite  straight-forward.  In  the 
second  case,  the  nonlinearity  is  considered  to  be 
"sandwiched"  between  the  two  elements  of  the  plant. 
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This  is  the  situation  found  most  often  in  physical 
systems.  All  normal  attempts  at  analysing  this 
system  fail  and  the  author  knows  of  no  text  or  ref¬ 
erence  which  has  developed  a  method  for  compensating 
this  type  of  a  system. 

1-2  THE  NONLINEARITY 

The  type  of  nonlinearity  which  will  be  con¬ 
sidered  is  shown  in  Fig,  1-2,  For  small  signals, 
the  operation  is  linear.  For  larger  inputs,  the 
system  operates  over  a  region  of  the  nonlinearity 
where  the  effective  gain  is  reduced.  In  this 
region  the  nonlinearity  is  assumed  to  be  a  straight 
line  with  slope  K,  where  K  is  less  than  1. 
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The  gradual  transition  between  the  two  segments 
was  included  to  cover  the  general  case  in  which 
the  transition  may  not  be  so  abrupt.  The  methods 
developed  in  this  thesis  are  completely  general 
in  that  the  gradual  transition  can  be  any  shape 
and  the  slope,  K,  can  be  varied  from  0  to  1  and 
the  process  is  still  valid. 

1-3  THE  DIGITAL  CONTROLLER 

The  digital  controller,  D(z),  is  the  compen¬ 
sating,  pulsed-data  network  which  would  normally 
be  used  with  a  linear  system  to  obtain  deadbeat 
response.  The  design  of  and  the  procedure  for 
simulating  the  controller  will  be  outlined  in 
the  following  chapter. 

1-4  THE  FORWARD -COMPENSATING  NETWORK 

The  forward-compensating  device,  Df(z),  is 
the  most  important  part,  and  the  object,  of  this 
thesis.  This  network,  which  in  general  contains 
nonlinear  elements,  senses  the  magnitude  of  the 
input  signal  and  adjusts  the  various  parameters 
of  the  digital  controller  so  that  deadbeat  response 
is  obtained  regardless  of  the  magnitude  of  the 
input  signal. 
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CHAPTER  2 


THE  LINEAR  SYSTEM 

In  the  design  of  the  digital  controller  for 
the  linear  system,  the  state  transition  approach 
shall  be  used.  Using  this  approach,  the  state 
transition  equations  are  used  to  describe  the 
state  of  each  variable  at  any  given  time.  This 
method  insures  that  there  will  be  no  inter-sampling 
ripple  and  that  the  response  is  indeed  deadbeat. 

The  closed-loop,  linear  system  is  shown  below 
in  Fig.  2-1. 


Fig.  2-1 


Sampling  period  =  T 
G(s)  =  1  /  s(s+l) 

D(z)  =  Kq  +  K1z"1+.  .  ,+KjjZ 
When  drawing  the  transition  flow  graph,  the 
digital  controller  can  be  considered  as  a  variable 
gain  network  as  discussed  in  Chapter  1»  The  zero- 
order  hold  measures  the  output  of  D(z)  at  the 
beginning  of  the  sampling  period  and  holds  this 
value  through  one  sampling  period.  The  transition 
flow  graph  is  shown  in  Fig.  2-2. 
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Fig,  2-2 

Using  Mason's  Gain  Formula,  the  state  transition 
equations  can  be  written.  By  taking  the  inverse 
Laplace  transform,  these  equations  can  be  written 
in  the  time  domain ,  This  set  of  independent 
equations  describes  the  response  of  the  state 
variables  during  a  certain  period  of  time.  For 
the  general  system,  having  n  state  variables,  n 
conditions  must  be  satisfied  by  n  equations.  The 
only  unknowns  in  the  system  are  KQ ,  K-^, 

Thus  n  unknowns  are  required  to  solve  the  n  simul¬ 
taneous  equations.  This  implies  that  n  sampling 
periods  are  required  to  obtain  a  specified  response 
from  an  'n'th  order  system.  Thus,  in  the  above 
example,  the  minimum  time  required  to  obtain 
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deadbeat  response  is  two  sampling  periods.  The 
two  equations  are  continuous  for  a  specified  'n', 
throughout  the  'n+1'  st  sampling  period.  Thus  the 
state  of  the  system  is  fully  described  by  the 
following  equations. 

x1(t)  =  (1  -  K(t0)(Jt-t0)  -  1  +  e" ^t:" to^ )  XiCto) 
+  (1  -  e  ^  to^)  X2 ( tQ) 

+  K(t0)  ( (t-tQ) -l+e~ ^ t“to^ )  r(t0) 

x2 ( t)  =  -K(t0)  (1  -  e  ^  to^) 

+  e"(t-to)x2(t0)  +  K(t0)(l-e"(t't°))  r(t0) 
For  tQ  =  kT,  t  =  (k+l)T,  r(t)  =  RQ  and  =  the 
gain  of  D(z)  during  the  'k+l'th  sampling  period; 

xi((k+l)T)  =  (1  -  Kk(T-14e'T))  xx(kT) 

(2-1)  K 

+  (l-e‘T)  x2(kT)  +  Kk(T-l+e'T)  RQ 

x2((k+l)T)  =  -Kk(l-e"T)  x^(kT)  +  e  Tx2(kT) 

(2-2)  _m 

+  Kk(l-e  X)  Rq 

The  conditions  for  deadbeat  response  to  a 
step  input  of  magnitude  R0  are: 

(a)  The  output,  x-^(t),  must  equal  RQ. 

(b)  All  other  state  variables  must  be  zero. 

It  is  required  that  the  system  approach  these 
conditions  within  as  few  sampling  periods  as 
possible.  Thus,  for  k  =  0,  with  all  initial 
conditions  being  zero,  it  is  required  that: 
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x1(T)  =  Kq (T  -  1  +  e"T)  R0  =  R0 
x2(T)  =  Kq  (  1  -  e'T)  R0  =  0 

From  these  equations,  it  is  verified  that  it  is 
not  possible  to  obtain  deadbeat  response  in  one 
sampling  period.  Assuming  that  deadbeat  response 
can  be  obtained  in  two  sampling  periods,  the  above 
values  of  x-^(T)  and  ^(T)  can  be  substituted  into 
equations  2-1  and  2-2  as  initial  conditions.  Then 
letting  k  =  1: 

Xl(2T)  =  (1  -  K1(T-l-te"T))K0(T-l-te‘T)  Rq 
+  (1  -  e”T)  Kq  (1  -  e'T)  Rq 
+  K.(T  -  1  +  e"T)  R0  =  RQ  and 
x2(2T)  =  -K^l  -  e”T)  Kq  (T  -  1  +  e”T)  R0 
+  e’T(l  -  e"T)  K0  Rq 
+  Ki  (X  -  e'T)  R0  =  0. 

Solving  these  equations  simultaneously  yields: 

K0  =  1  /  T(l-e“T) 
and  Kj^  =  1  /(l+T-e+T) 

If  a  sampling  period  of  one  second  is  used,  the 
values  of  KQ  and  are: 

Kq  =  1.5820 
and  Kx  =-1.3922 
xx(T)  =  R0/  (e-1)  =  .5820  R0 

e  (T)  =  R0  -  x^T)  =  .4180  RQ 

h(T+)  =  Kx  e (T)  =  .4180  (-1.3922)R0  =-.5820  R0 
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Thus  : 


D(z)  = 


1.5820  -  .5820  z"1 

1  +  .4180  z"1 


.2-1  SIMULATION 

Once  D(z)  has  been  calculated,  the  system  can 
be  simulated  on  the  analog  computer.  To  simulate  the 
digital  controller,  the  circuit  shown  in  Fig.  2-3  is 
most  practical.  The  design  and  operation  of  the 
controller  is  described  in  Appendix  A.  With  reference 
to  Fig.  2-3,  the  initial  error  signal,  e(0+),  which  is 
equal  to  RQ ,  is  multiplied  by  K0  and  held  by  the  zero- 
order  hold  for  one  sampling  period.  At  the  end  of  one 
sampling  period,  the  output  of  the  plant,  x^(t),  will 
have  reached  a  value  of  .5820  RQ.  Thus  e(T)  is  equal 
to  .4180  R0.  Since  the  potentiometer,  P2 ,  is  set  at 
a  value  of  .4180,  then  the  output  of  amplifier  A  will 
be  zero  at  the  beginning  of  the  second  sampling  period. 
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Fig.  2-3 
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Since  the  potentiometer,  ,  is  set  at  .5280,  the  out¬ 
put  of  amplifier  B  will  be  equal  to  -.5280  RQ  at  the 
instant  T+.  The  zero-order  hold  then  holds  this  value 
for  one  sampling  period  and  the  system  responds  to 
this  negative  step.  At  the  end  of  the  second  sampling 
period,  the  output  of  the  system  will  be  equal  to  RQ 
and  the  state  variable,  X2(t),  will  be  equal  to  zero. 
Thus  deadbeat  response  has  been  achieved.  It  should 
be  noted  that: 

(a)  since  the  output  of  Amplifier  A  is  zero  at 
the  instant,  T+,  the  output  from  the  time  delay  will 
be  zero  for  all  time  after  2T. 

(b)  since  the  output  of  the  plant  is  equal  to 
R  at  t  =  2T,  the  error  signal  will  be  zero  for  all 
time  after  2T. 

The  complete  computer  schematic  for  this  linear 
system  is  shown  in  Fig.  2-4. 

2-2  RESULTS 

Using  the  network  shown  in  Fig.  2-4,  it  is 
desirable  to  observe  whether  the  system  response 
meets  the  specifications  required  for  deadbeat  response. 
In  order  to  compare  the  results  which  will  be  obtained 
for  the  compensated  systems,  it  is  also  necessary  to 
observe  some  measure  of  sensitivity  for  this  linear 
system.  To  do  this,  certain  parameters  can  be  adjusted. 
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QOQ  =  Ro/100 

POl  =  „  4180 

P04  =  K  /2  =  „  7910 

Flg°  2"4  Q04  =  . 5820 

The  percentage  deviation  from  the  desired  response  can 

then  be  measured.,  If,  for  example ,  in  the  actual  plant, 
the  sampling  period,  T,  and  the  time  constant, T  s  are 
not  known  to  any  degree  of  accuracy,  then  the  digital 
controller  will  not  yield  perfect  deadbeat  response. 

The  percentage  overshoot  (or  undershoot)  for  various 
changes  in  T  and  TT  are  shown  in  graphs  2-1  and  2-2. 

The  results  indicate  that  if  the  sampling  period 
is  varied  by  10%  from  the  value  for  which  it  was  designed, 
the  response  would  be  in  error  by  no  more  than  97o.  There 
is  an  almost  linear  relationship  between  the  percentage 
change  in  the  sampling  period  and  the  percentage  over¬ 
shoot  of  the  response. 
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Similar  results  occur  for  variations  in  the  time 
constants  X*  The  sensitivity  is  almost  constant  at  0.5, 
meaning  that  for  a  percentage  change  of  10%  in  the  time 
constant,  the  percentage  overshoot  is  approximately  57,. 

2-3  CONCLUSIONS 

The  results  illustrate  ho w  the  digital  controller, 
used  in  the  linear  system,  can  be  designed  to  obtain  minimal 
response  to  an  accuracy  of  less  than  1%.  Although  this 
system  is  not  insensitive  to  parameter  variations,  a 
certain  amount  of  variation  can  be  allowed  and  the 
deterioration  in  the  response  is  not  drastic.  Thus  if, 
for  example,  the  designer  is  willing  to  tolerate  a  107, 
overshoot  or  undershoot  in  the  response,  then  the  samp¬ 
ling  period  can  be  allowed  to  drift  by  as  much  as  107o. 

In  practical  applications,  the  sampling  period  can 
usually  be  kept  accurate  to  within  0.1  to  0.01  percent. 
Similarily,  the  time  constant  of  the  system  needs  to  be 
known  only  to  within  approximately  107,  of  its  actual 
value  in  order  to  design  a  controller  which  will  yield 
deadbeat  response  to  within  5%  accuracy. 
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CHAPTER  3 


THE  NONLINEAR  SYSTEM 
WITH  ISOLATED  NONLINEARITY 

The  nonlinearity  which  was  considered  cannot,  in 
general,  be  described  exactly  by  any  mathematical 
expression.  The  analysis  of  the  system  response  for 
various  magnitudes  of  inputs  cannot  be  done  by  using 
strictly  mathematical  approaches.  An  attempt  was  made 
to  approximate  the  curve  (shown  in  Fig.  3-1)  using  an 
expression  of  the  form: 

eQ  ®  (ei  +  aKe?)  /  (1  +  ae.)  , 
and  use  a  complimentary  function  in  the  compensating 
network.  This  approach  was  abandoned  for  the  following 
reasons . 

(a)  The  difference  between  the  actual  curve  and 
the  curve  described  by  the  approximate  function  could, 
in  some  cases,  make  the  use  of  this  method  impractical. 
Because  of  this  there  is  a  loss  in  generality. 

(b)  If  the  actual  curve  did  indeed  correspond 
exactly  to  the  approximation,  the  error  introduced  in 
compensating  the  system  with  the  complimentary  math¬ 
ematical  expression  was  in  some  cases  too  large. 

3-1  COMPENSATION 

To  compensate  for  the  nonlinearity,  consider  the 


system  shown  in  Fig.  3-2, 
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Fig.  3-2 

In  the  previous  chapter:  it  w as  shown  that  the 
input  to  the  plant ,  h(t),  is  composed  of  two  pulses, 
each  of  duration  T  seconds.  The  magnitude  of  each 
pulse  is  linearly  dependent  upon  the  amplitude  of  the 
input  signal,  R0 .  The  nonlinearity,  in  this  case, 
affects  the  magnitude  of  the  output  of  the  digital 
controller.  The  input  signal  to  the  nonlinearity 

i 

must  be  increased  to  such  a  value  that  the  output,  h(t), 
is  the  required  magnitude.  This  is  done  by  multiplying 

KoRo  by  Ko  and  eiKiR0  by  Ki‘ 

For  the  first  sampling  period: 

ho  -  KoRo  =  f<ho> 
h*  =  KoR0k|  =  h0K^ 

Ko  =  h*/ho  =  §o(Ro) 

Similarily  for  the  second  sampling  period: 
h1  =  R0K1e^  =  f(h^) 
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h*  =  KieiR0Kf  =  hj_Kf 
K1  =  hl  1  hl  =  SfRo> 

£  f 

Kq  and  can  then  be  plotted  as  a  function  of  RQ 
The  results  are  shown  on  graph  3-2. 

3-2  SIMULATION 

Using  the  results  of  Sec.  3-1,  it  is  possible  to 
simulate  a  system  to  test  the  validity  of  this  feed¬ 
forward  method  of  compensation.  This  system  is  shown 
in  Fig.  3-3.  A  similar  system,  with  no  feedforward 
compensation,  was  simulated  to  illustrate  the  effect¬ 
iveness  of  the  compensation.  The  results  are  shown 
in  Graph  3-1. 


Fig.  3-3 
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Since  the  two  graphs,  and  K][  versus  R0,  are 

identical,  except  for  the  scale  factor  on  the  horizontal 

axis,  one  Diode  Function  Generator  (D.F.G.)  can  be  used 

to  generate  both  functions.  The  sensitivity  of  this 

system  can  then  be  measured  using  the  same  method  as 

see 

in  Chapter  2,  to^if  there  has  been  any  degrading  of 
the  overall  response  due  to  the  addition  of  the  compen¬ 
sating  network.  The  results  are  shown  on  graphs  3-3 
and  3-4. 

3-3  APPROXIMATE  METHOD  OF  COMPENSATION 

In  practical  applications  of  this  method,  it  is 
often  impractical  and  very  expensive  to  use  D.F.G.'s 
to  generate  the  nonlinearity  used  in  the  feedforward 
network.  This  problem  can  be  overcome  by  approximating 
the  nonlinearity  with  two  straight  line  segments.  This 
function  is  quite  simple  to  generate  and  the  degradation 
of  the  response  is  not  drastic.  The  network  used  to 
generate  the  nonlinear  function  is  shown  in  Fig.  3-4. 

If  this  function  is  substituted  for  the  nonlinear 
element  in  the  feedforward  network  in  Fig.  3-3  (page  19), 
that  network  can  be  used  to  measure  the  response  for 
various  magnitudes  of  step  inputs  and  to  determine  the 
degree  to  which  the  response  has  deteriorated  due  to 
the  use  of  this  approximate  method.  The  results  are 
shown  in  Grapji  3-5. 
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P08  =  0. 1000 
Q08  =  0.  1122 
P09  =  0.  38  20 


The  accuracy  of  the  response  has  not  been  reduced 

drastically.  In  the  range,  7<  RQ<  11,  the  function, 

f  f 

Kq  vs.  Rq  ,  is  in  error  by  the  greatest  amount.  KQ  is 

less  than  the  required  value  and  as  a  result,  hQ  is 

too  small  and  the  response  is  an  undershoot. 

3-4  CONCLUSIONS 

The  results  of  this  chapter  illustrate  the  necess¬ 
ity  for  compensation  and  also  ho w  a  system  with  an 
isolated  nonlinearity  such  as  the  one  used  in  this 
example,  can  be  compensated  so  that  almost  perfect 
deadbeat  response  is  obtained,  regardless  of  the  mag¬ 
nitude  of  the  step  input.  The  effect  of  changes  in 
the  sampling  period,  T,  and  the  time  constant,TT  ,  is 
no  more  serious  than  for  the  linear  case. 

It  should  be  noted  that  the  nonlinearity  used 
in  this  case  is  one  of  the  worst  possible  situations 
one  could  expect. 
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The  saturation  type  of  nonlinearity  used  here  has  unity 
gain  in  the  linear  region  and  in  the  saturation  region, 
the  slope  of  the  nonlinearity  is  0.1.  If  this  slope 
was  zero,  then  this  procedure  could  not  be  considered 
for  compensation.  Any  other  methods  which  would 
require  only  two  sampling  periods  would  also  fail. 

One  must  then  resort  to  using  three  sampling  periods 
to  obtain  minimal  response. 

This  chapter  also  illustrates  how  an  approximate 
method  of  compensation  can  be  used  to  obtain  satisfac¬ 
tory  results.  Various  methods  of  approximating  the 
nonlinearity  can  be  used,  depending  upon  the  range 
over  which  the  system  is  to  operate  and  the  type  of 
error  which  can  be  tolerated. 

This  method  can  be  applied,  with  minor  changes, 
to  systems  in  which  the  nonlinearity  is  nonsymmetrical . 

By  using  one  unit- time-delay ,  it  can  also  be  adapted 
to  compensate  systems  where  successive  inputs  are  applied. 
This  includes  either  positive  or  negative  inputs. 
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CHAPTER  4 


NONLINEAR  SYSTEM  WITH  INTEGRATED  NONLINEARITY 
In  this  chapter,  as  in  Chapter  3,  a  simple,  type  1, 

second  order  system  will  be  considered.  The  nonlinearity 

will  be  considered  to  be  an  integral  part  of  the  plant. 

most  texts  or  references  on  the  subject,  the  nonlinearity  i§ 

considered  to  be  separable  from  the  remainder  of  the 
plant.  In  many  practical  situations  however,  this  is 
not  so.  The  block  diagram  of  the  system  to  be  consid¬ 
ered  is  shown  in  Fig.  4-1. 


Fig.  4-1 

For  this  system,  the  overall  transfer  function 
can  no  longer  be  described  by  simple  linear  equations. 
The  input-output  relation  for  each  element  must  be 
considered  separately. 
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(a)  h(t)  is  a  series  of  steps  formed  by  using 
the  zero-order  hold  to  hold  the  values,  h(01'),  h(T+) , 
h(2T+)  ,  .  ..  s  h(nT+)  through  the  'n+l'st  sampling 
period . 

(b)  ^(t)  is  the  exponential  response  to  the 
step  function,  h(t). 

(c)  x^(t)  is  the  output  of  the  nonlinearity. 

This  response  is  a  .'distorted  exponential'  function, 
the  relation  between  X2(t)  and  X2(t)  being  determined 
by  the  nonlinearity. 

(d)  The  output  of  the  system,  x-^(t),  is  the 
integral  of  X2(t)  or  simply  the  area  under  the  curve, 
X2(t)  vs.  time. 

4-1  COMPENSATION 

It  is  obvious  that  the  plant  response  is  dependent 
upon  Rq,  the  magnitude  of  the  step  input,  and  as  a 
result,  the  operation  of  the  Digital  Controller  must 
be  analysed.  In  Chapter  2,  it  was  shown  that  the 
transfer  function  of  a  digital  controller  can  be  des¬ 
cribed  by  the  following  equation. 

D(z)  .  h(0+)  +  h(T+) 

I  +  e(T+) 

In  this  system  there  are  three  controlling  factors, 
h(0+) ,  h(T+)  and  e(T+).  For  large  values  of  RQ,  the 
two  values,  h(0+)  and  h(T+),  will  have  to  be  multiplied 
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by  certain  constants  to  make  up  for  the  loss  in  gain 
introduced  by  the  nonlinearity.  It  is  also  reasonable 
to  assume  that  the  output,  x^(T),  will  not  be  equal  to 
0.5820  R0  and  as  a  result,  e(T+)  will  not  equal  .418R0. 
This  parameter  will  also  have  to  be  adjusted  to  diff¬ 
erent  values  for  different  magnitudes  of  inputs.  Thus 
it  appears  as  though  the  three  parameters  of  the 
Digital  Controller  will  all  be  nonlinear  functions  of 
R  ,  the  magnitude  of  the  step  input.  To  determine  a 
method  to  calculate  these  nonlinear  functions,  the 
operation  of  the  system  must  be  considered  in  a  step- 
by-step  fashion.  Certain  restrictions  will  then  be 
applied  and  the  unknown  values  can  then  be  obtained. 

Referring  to  Fig.  4-1,  the  typical  response  curves 
of  the  various  element  outputs  can  be  sketched  (Fig.  4-3). 
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Normally,  the  relationship  between  X2(t)  and  X2(t) 
cannot  be  described  by  a  simple  mathematical  expression. 
However,  if  the  nonlinearity  is  made  up  of  (or  approx¬ 
imated  by)  straight  line  segments,  then  separate  equa¬ 
tions  can  be  written  to  describe  this  relationship 
through  each  section  of  the  nonlinearity.  The  sets  of 
points  used  in  making  up  the  nonlinear  functions,  to 
be  used  in  the  compensating  network,  can  be  calculated 
by  repeating  the  following  step-by-step  procedure. 


Fig.  4-4 
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(a)  Assume  that  h(0+)  is  of  such  magnitude  that 
X£(T)  is  equal  to  An.  (Refer  to  Fig.  4-4) 

Then  A^  =  h(0+)  (1  -  e‘T)  and  h(0+)  =  A^/ (1  -  e'T) 

(b)  X2(2T)  =  0  is  one  condition  which  must  be 
satisfied  in  order  to  have  deadbeat  response. 

Thus:  h(T+)  =  An  -  h(0+) 

(c)  Calculate  the  times  required  for  X£(t)  to 
rise  to  each  of  the  values  A2>  A3,  ..o,  An.  That  is: 

h(0+)  (1  -  e'tl)  =  A2 
e'6!  =(h(0+)  -  A2)/h(0+) 
t,  =  -antilog(  l-A2/h(0+)) 

(d)  Similarly  s  calculate  the  time  required  for 
X2(t)  to  decay  to  An_^J  ^-2 5  •••> 

(e)  Determine  the  equation  for  the  straight  line 
segments  from  A^B^  to  A2B2,  A2B2  to  A3B3,  ...,  and 
from  An_^Bn_-^  to  AnBn.  This  gives  the  relationship 
between  x^(t)  and  X2(t)  for  the  time  intervals 
calculated  previously. 

(f)  Using  these  equations,  integrate  X2(t) 

through  each  of  the  time  intervals.  The  sum  of  these 

/2T  * 

x0(t)  dt  which  is  equal  to 
u 

x^(2T)  which  is  required  to  be  equal  to  R0. 

/T  * 

X2(t)  dt  which  equals  x^(T). 

From  this,  e(T)  =  RQ  -  x^(T). 
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(h)  R0K0K^  =  h(0+)  =  h(0+)/  R0Kq 

R0e1K1Kf  =  h(T+)  k£  =  h(T+)/  R^^ 

RQe^  =  e(T+)  el(R0>  =  e(T+)/  R0 

These  values  form  one  set  of  points  for  the  three 
curves  relating  K^,  and  e^  to  RQ,  and  are  valid  for 
only  one  value  of  R0,  This  whole  series  of  calculations 
can  be  done  very  simply  by  using  the  digital  computer. 

The  computer  program,  along  with  the  explanation  of  the 
notations  used,  is  shown  in  Appendix  B0 

x: 

The  results  of  this  program  give  the  values  of  K£, 
k|  and  e^  required  for  various  values  of  R0.  The  values 
of  and  k£  were  found  to  be  equal  for  each  value  of  RQ, 
This  means  that  the  outputs  of  the  original  digital  con¬ 
troller,  and  e^K^RQ,  need  to  be  multiplied  by  the 

same  correction  factor,  K^,  for  each  value  of  RQ.  The 
functions,  and  e^  vs.  RQ,  are  plotted  on  Graph  4-1. 

4-2  SIMULATION 

Using  the  above  results,  a  system  can  be  simulated 
to  test  the  response  for  various  magnitudes  of  step 
inputs.  The  system  is  shown  in  Fig.  4-5. 

4-3  RESULTS 

Using  the  system  shown  in  Fig.  4-5,  the  response 
to  various  magnitudes  of  step  inputs  was  measured.  To 
illustrate  the  effectiveness  of  this  feedforward  method 
of  compensation,  the  same  system  with  no  compensation 
was  tested. 
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Fig„  4"5 

Th,e  nonlinear  functions  used  in  the  compensating 
network  were  then  approximated  using  an  arrangement 
similar  to  the  one  shown  in  Fig.  3-4.  The  results 
of  these  three  tests  are  shown  in  Graph  4-2. 

4-4  CONCLUSIONS 

Results  of  this  chapter  demonstrate  how  this 
feedforward  method  of  compensation  can  be  used  to 
compensate  for  a  nonlinearity  which  cannot  be  isolated 
from  the  remainder  of  the  plant.  The  results  of  the 
compensated  system  compare  favorably  with  those  of 
the  linear  system. 


-  ..  - — •  —  -  * 


oox\„j 


! 


■'1 


'• 


* 


i~±  .ail 


‘  .  .  -  S  ' 

§n  ^fii  nscftioo  sift  nl  beau  auoisamj^  ix-^rrilr  yn  sdT 
. 

e-l'  Id  3.01  arfT  ,£-£  ♦  gi’3  tfi  riworia  ono  arts  o3  iBltmi  a 

-)  ni  r 

HT1Gt£JJOViO) 

ax/to  wori  noireb  ’Ja3q3rii>'€-tii3  lo  2J Ju  etefl 

. 


.  ..  ioai  *  3  •  ■■■•-.:  >  v^l'r&s.iHno . ft  vox  9^sftiz9,qmoo 

#  • 

sri3  . Snsiq  adj  io  labni Bros’*  srta  meal 


. 


'4 


Ul.  . '  •  -  ;  :  •:) 

m*  ve  so*  r  toHt^ 
rtrwRj  irioH  smutest;  : 


■  ■■ 

■?>? 


\j 


■ 

Z1  *  -  :  ■ 


> 


rimw  wataye 


WOITAe-ViSiiCD  OM 


37 


The  response  obtained  using  no  compensation  is 
sufficiently  accurate  for  the  smaller  values  of  R0. 
However,  as  the  magnitude  of  the  step  input  increases, 
the  response  continually  deteriorates  until  it  does  not 
even  approximate  minimal  response. 

Although  the  results  obtained  using  the  straight 
line  approximations  are  degraded  slightly  in  comparison 
to  the  correctly  compensated  system,  they  illustrate 
how,  in  practical  situations,  this  approximate  method 
of  compensation  is  quite  satisfactory  if  one  is  willing 
to  accept  'almost  perfect'  deadbeat  response. 

With  a  few  minor  changes,  the  system  shown  in 
Fig.  4-5  can  be  adapted  to  compensate  for  the  case 
where  successive  step  inputs  are  applied.  This  method 
of  compensation  can  also  be  applied  to  the  case  where 
negative  steps  are  applied  and  the  cases  where  the 
nonlinearity  is  either  symmetric  or  nonsymmetric . 

These  adaptations  will  be  outlined  more  fully  in 
Chapter  7. 
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CHAPTER  5 


DESIGN  OF  COMPENSATING  NETWORKS 
USING  AN  EXPERIMENTAL  METHOD 
The  results  of  Chapter  4  illustrate  how  effective 
the  feedforward  method  is  in  compensating  for  certain 
types  of  nonlinearities  in  sampled-data  systems.  The 
only  problem  arises  in  deriving  the  nonlinear  functions 
which  are  to  be  used  in  the  feedforward  network.  In 
most  cases ,  a  digital  computer  may  not  be  available  to 
do  the  extensive  calculations  required.  If  the  input- 
output  characteristic  of  the  plant  nonlinearity  is  not 
known  or  if  the  time  constant  in  the  l/(s+'C)  term  can¬ 
not  be  determined  exactly,  then  the  method  outlined  in 
Chapter  4  cannot  be  used. 

An  alternative  method  of  deriving  the  'feedforward 
functions’  will  be  outlined  in  this  chapter.  The 
experimental  method  used  here  is  simplified  by  utilizing 
the  'Repetitive  Operation'  (Rep.  Op.)  Mode  of  the 
analog  computer.  High  speed  electronic  comparators 
are  used  to  apply  a  series  of  steps  to  the  open  loop 
system.  The  basic  circuit  used  is  shown  in  Fig.  5-1. 
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5-1  REPETITIVE  OPERATION  MODE 

The  Rep.  Op.  mode  of  the  analog  computer  is  a  special 
mode,  in  which  the  response  of  any  circuit  is  speeded  up 
by  a  factor  of  500.  The  computer  automatically  switches 
from  ’initial  Conditions’  Mode  to  ’Operate’  Mode.  The 
operate  time  can  be  selected  as  10,  20,  40,  or  80ms. 
duration . 

5-2  TIME  SCALING 

In  this  example,  the  system  which  is  simulated  on 
the  computer  can  be  slowed  down  by  a  factor  of  20,  so 
that  in  Rep.  Op.,  one  second  of  real  time  corresponds 
to  20/500  or  40ms.  of  computer  time.  The  entire  response 
of  the  system,  corresponding  to  two  sampling  periods, 
can  then  be  displayed  in  80  ms.  Including  reset  time, 
the  computer  takes  approximately  100  ms.  of  real  time 
for  each  response.  The  response  is  then  repeated  at  a 
rate  of  about  10  times  per  second.  The  various  outputs 
can  then  be  displayed  on  a  C.R.T.  and  appear  as  continuous 
curves . 

5 “3  RESULTS 

Using  the  circuit  shown  in  Fig.  5-1,  with  the  com¬ 
puter  in  Rep.  Op,,  and  the  various  outputs  displayed  on 
the  C.R.T. ,  the  two  potentiometers,  P01  and  P02 ,  can 
then  be  adjusted  until  x^Ct)  exhibits  deadbeat  response. 
This  means  that  X2(2T)  must  equal  zero.  The  magnitudes 
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of  the  pulses,  h(0+)  and  h(T+),  as  well  as  x-^(T)  can 
then  be  measured.  This  process  is  then  repeated  for 
several  values  of  X2(2T).  With  the  values  of  hQ ,  h^, 
x^(T)  and  x-^(2T)  recorded  for  each  case,  Kq,  and 
e(T)  can  then  be  calculated. 

h(0+)  =  KqR0k|  k|  =  h(0+)/KoRo 

h(T+)  =  e1K1R0Kf  =  h(T+)/e1K1R0 

e(T+)  =  eiR0  e  x (R0 )  =  e(T+)/R0 

The  results  of  Chapter  4  indicate  that  Kq  should 
be  equal  to  K^.  The  average  value  can  be  used  in  this 
case  to  obtain  a  value  of  K^.  The  actual  computer  dia¬ 
gram  used  to  illustrate  this  experimental  method  is 
shown  in  Fig  5~2,  the  results  being  displayed  on  Graph  5-1. 
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To  determine  how  valuable  this  method  is,  the  response 
of  the  system  was  measured  using  the  experimentally-obtained 
functions  in  the  feedforward  loop.  The  response  is  shown 
in  Graph  5-2  and  compared  with  the  response  obtained  using 
the  theoretical  functions  in  the  feedforward  loop  of  the 
system. 

5 “4  CONCLUSIONS 

The  overall  response  of  the  system,  using  the  exper¬ 
imentally  obtained  results,  has  deteriorated  by  approx¬ 
imately  a  factor  of  two.  Much  of  this  error  can  be  reduced 
in  the  linear  region,  since  the  designer  would  normally 
know  the  characteristic  of  the  nonlinearity  in  the  linear 
region  and  would  use  a  value  of  one  for  Kj .  If  the  designer 
intends  to  use  straight  line  approximations  to  generate  the 
nonlinear  functions,  then  these  experimentally-obtained 
results  can  be  used  as  a  guide  in  generating  the  functions. 

Although  the  method  outlined  in  this  chapter  can  be 
applied  to  the  system  studied  in  Chapter  3,  it  would  be  of 
very  little  value  since  that  was  a  relatively  simple  system. 
Even  in  cases  where  a  digital  computer  is  available,  this 
method  is  very  effective  since  in  general,  the  nonlinearity 
is  not  known  to  a  high  degree  of  accuracy  and  the  variation 
of  the  characteristic  is  usually  such  that  the  response 
obtained  using  the  experimentally-obtained  results  is  just 
as  accurate  as  that  obtained  using  the  theoretical  results. 
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CHAPTER  6 


OPTIMIZATION  OF  THE  SYSTEM  USING 
THREE  SAMPLING  PERIODS 

In  the  systems  which  have  been  studied  in  the 
previous  chapters ,  the  range  of  inputs,  RQ ,  was 
appropriately  chosen  so  that  the  input  to  the  non¬ 
linear  element  covered  a  convenient  range  of  the  non¬ 
linear  region.  In  practical  systems,  the  nonlinearity 
may  be  such  that  the  saturation  becomes  very  severe  and 
the  output  may  even  approach  a  finite  maximum  (see  Fig.  6-1). 
This  situation  would  require  excessively  large  plant  inputs 
in  order  to  obtain  the  required  output.  In  cases  of  complete 
saturation,  the  required  output  may  never  be  obtained. 


period.  This  however  degrades  the  response  for  smaller 
magnitudes  of  inputs.  The  other  obvious  solution  is  to  use 
three  sampling  periods  for  the  larger  inputs  and  the  regular 
two  periods  for  the  smaller  inputs. 
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Since  using  three  sampling  periods  necessarily  means 
the  addition  of  one  variable,  a  number  of  alternatives  are 
available  in  the  design  of  the  digital  controller.  In 
general,  one  would  normally  try  to  optimize  one  or  more 
items  in  the  system  response,  the  item  to  be  optimized 
depending  upon  the  system  under  consideration. 

6-1  SOLUTION 

In  this  chapter,  the  system  used  in  Chapter  4  will 
again  be  considered  (see  Fig.  6-2). 


Fig.  6-2 


If  the  power  input  to  the  plant  is  proportional  to 

2 

h  (t),  then  one  would  normally  try  to  maximize  the  output, 

R0 ,  with  respect  to  the  input  power.  Since  the  output  is 

* 

equal  to  the  integral  of  the  state  variable,  X2(t),  and 
since  X2(t)  is  related  to  X2(t)  through  the  nonlinear 
element,  consider  first' the  integral  of  X2(t)  for  various 
values  of  h(t)  (see  Fig.  6-3).  In  the  range  0  <  t  <  T: 

x2(t)  =  hQ(l  -  e  fc) 

A1  =  h0/  C1  "  e” t)dt  =  hQ(t  +  e_t)]^  =  h0 (T  +  e  T-  1) 
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In  the  range  T  <  t  <  2T: 

x2(t)  =  x2(T)  +  (h-j^  x2(T))(1  -  e"  ( t-T) ) 

=  hx(l  -  e'^-T))  +  hQ(l  -  e“T)e-^-T> 
x2(2T)  =  (hx  +  hoe"T)(l  -  e"T) 

A2  ~  x2^fc)  dt  =  dl(T  +  e_T“  1)  +  H0(l  "  e-T)2 

T  y 

In  the  range  2T  <  t  <  3T: 

x2(t)  =  x2(2T)  +  (h2-  x2(2T))(1  -  e  ^  2T^) 

~  h2(l  -  e-(t-2T))  +  (h-^4-  hQe  T)(l  -  e  •r)e"^t"^T^ 
However ,  the  requirement  for  deadbeat  response  is  that 
x2(3T)  =  0.  Thus: 

h2(l  -  e"T)  +  (hx  +  hoe"T)e"T(l  -  e"T)  =  0 

h2  =  -(h^  +  hoe-T)  e  T 

x2(t)  =  (h-L  +  hoe“T)  (  e~  (t_2T)  _  e  T) 

A3  =  2rf x2 ^ dt  =  +  h0e"T)  (!  -  Te  T-  e  T) 

The  total  power  is  Arj,  =  A-^  +  A2  +  A^ 

=  h0T(l  -  e"2T)  +  hxT(l  -  e"T)  Eqn .  6-1 

The  total  power  into  the  plant  is  assumed  to  be  proportional 
to  (h^  +  h2  +  h|)  =  h^  +  h|  +  (hi  +  hoe”^)^e  2^  oi  P  Eqn.  6-2 
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Using  equations  1  and  2,  the  relation  between  h-^  and  hQ  can 
be  found  so  that  P  is  a  minimum  for  a  constant  area,  A^. 
From  Eqn.  6-1: 

•-  C  -  hiT(l  -  e"T) 
ho  -  f<hl>  -  - T  a  .  e-2f) 

dh0  =  .  (  1  -  e~T) 

dhj_  (  1  -  e'2T) 

From  Eqn .  6-2 : 


=  0 

dhi  ()hi  dh^  ^hQ  dh^ 

Solving,  h-^  =  h  (  1  -  e  ) 

This  is  the  condition  to  be  satisfied  in  order  that,  for  a 
specified  value  of  A^,  the  input'  power  is  minimized.  Alter¬ 
natively,  for  a  specified  power  input,  the  maximum  output 
is  obtained  if  h-^  =  hQ(l  -  e”^)  .  This  means  that  ^(t)  is 
constant  through  the  second  sampling  period.  That  is, 

X2 ( t)  =  hQ(l  -  e  for  T  <  t  <  2T. 

Typical  response  curves  are  shown  in  Fig.  6-4. 
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study  the  relationship  between  X2(t)  and  x^Ct) 
3y  3-p 

let  q  yx2(t)  dt  =  A  and  dt  =  A  .  For  a  given 


k 


value  of  A,  the  maximum  value  of  A  occurs  when  the  res- 

_  tp 

ponse  is  as  shown  in  Fig.  6-4  (i.e.  h^  =  hQ(l  -  e  )). 

This  means  that  x^(3T)  =  RQ  is  maximized.  In  other  words, 
for  a  specific  value  of  RQ ,  the  input  power  to  the  plant 

and  the  maximum  value  of  X2(t)  (and  also  x^(t))  can  both 

-  T 

be  minimized  by  applying  the  condition,  h^  =  hQ(l  -  e  ). 
6-2  SIMULATION 

Using  the  conditions, 
h-,  =  h  (1  -  e~^)  and 


o 


-Ts  -T 


h2  =  “^hl  +  hoe  )  e  =  "hoe 


-T 


hQ  +  hQ(l  -  e~ 1)z~ x  +  ( -hQe~T)  z~2 


D(z)  can  readily  be  written  as, 

»n  t  «0Vi-  -  *  Tn-~1 
1  +  e(T+)  z'1  +  e(2T+)  z‘2 

By  referring  to  Fig.  6-4a,  it  is  noticed  that  for 
the  system  using  three  sampling  periods,  the  area  under 
the  curve  of  X2(t)  vs.  time  is  exactly  double  the  area 
if  two  sampling  periods  were  used  (for  equal  magnitudes 
of  h(0+)).  Thus,  the  gain,  KQ,  required  using  three 
sampling  periods  is  %  of  that  required  when  using  two 
sampling  periods.  Thus: 

,  »  .7910  (  1.00  +  0.6320  z"1  -  0.3680  z"2) 

D(z)  =  - — — — ~ t - - - 7— 

(1.00  +  0.7090  z"1  +  0.2140  z~z) 

Using  the  above  design,  the  system  will  exhibit 


deadbeat  response  in  three  sampling  periods,  as  long  as 
the  operation  of  the  system  is  in  the  linear  region. 
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r(t). 


It  can  also  be  seen  from  Fig.  6-4(a)  that  the  operation 
of  the  system  is  linear  for  all  magnitudes  of  step  inputs 
from  0  to  20  volts.  For  inputs  larger  than  this,  feed¬ 
forward  compensation  as  discussed  in  Chapter  4  is  necessary 
If  the  operation  of  the  system  is  limited  to  inputs 
in  the  range,  0  to  20  volts,  then, the  following  system, 
shown  in  Fig.  6-5,  can  be  used. 
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6-3  CONCLUSIONS 

In  this  chapter  a  different  method  of  feedforward 
compensation  was  studied.  The  feedforward  channel  was 
used  only  to  switch  comparators  which  selected  one  of 
two  digital  controllers.  The  two  digital  controllers 
were  designed  so  that  the  response  of  the  system  was 
minimal  as  long  as  the  system  operated  over  the  linear 
portion  of  the  nonlinear  element.  One  digital  controller 
was  designed  so  that  deadbeat  response  was  obtained  in 
two  sampling  periods  for  any  step  input  from  0  to  10  volts. 
The  other  was  designed  for  using  three  sampling  periods 
for  any  inputs  from  0  to  20  volts. 

For  inputs  larger  than  20  volts,  one  could  use 
compensation  similar  to  that  discussed  in  Chapter  4. 
Depending  upon  the  range  of  inputs  used,  the  type  of  non¬ 
linearity  used  and  the  quality  of  response  required,  one 
may  use  any  combination  of  the  three  methods  mentioned. 

Unlike  the  systems  studied  in  Chapters  3  and  4,  the 
system  shown  in  Fig.  6-5  can  accept  any  number  of  succes¬ 
sive  step  inputs  provided  they  are  no  larger  than  20  volts 
and  they  are  not  applied  within  3  (or  two)  sampling 
periods  of  each  other. 
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CHAPTER  7 


ADDITIONAL  MODIFICATIONS  USING 
FEEDFORWARD  TECHNIQUES 

In  previous  chapters,  the  systems  which  were  studied 
were  considered  to  be  initially  at  rest.  The  compensation 
was  then  studied,  assuming  that  any  step  input,  over  a 
specific  range  of  magnitudes,  could  be  applied.  The  plant 
nonlinearity  was  also  assumed  to  be  symmetrical.  That  is, 
the  gain  characteristic  was  assumed  to  be  identical  for 
positive  and  negative  inputs.  No  consideration  was  given 
as  to  what  would  happen  if  the  input  was  suddenly  removed 
or  if  successive  inputs  were  applied  or  if  negative  inputs 
were  applied.  A  linear  system,  using  the  correct  digital 
controller,  will  respond  in  any  of  these  situations  with 
deadbeat  response,  A  nonlinear  system,  however,  in  many 
of  these  situations  will  not. 

In  this  chapter,  two  different  nonlinear  systems  will 
be  compensated  so  that  the  response  in  any  of  these 
situations  will  be  deadbeat. 

7-1  SYSTEM  WITH  ISOLATED  NONLINEARITY 
7-1-A  NONSYMMETRICAL  NONLINEARITY 

With  reference  to  Chapter  3,  the  system  shown  in 
Fig.  7-1  will  be  considered.  In  this  system,  the  non¬ 
linearity  was  symmetrical  and  as  a  result,  one  D.F.G. 
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could  be  used  to  generate  and  Kj .  With  a  nonsymme trical 
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Fig.  7-1 

nonlinearity ,  two  separate  nonlinear  functions  must  be 
generated ,  one  determined  by  the  positive  portion  of  the 
nonlinearity  and  the  other  depending  upon  the  negative 
portion.  The  feedforward  portion oof  the  system  would 
need  to  be  modified  as  shown  in  Fig.  7-2.  The  magnitude 
scaling  has  been  eliminated  for  simplification.  D.F.G. 

'A!  is  the  nonlinear  function  generator  which  generates 
the  factor  by  which  the  first  positive  step,  h(0+),  is 
multiplied  (assuming  a  positive  step  input  is  applied 
to  the  system).  Its  characteristic  is  determined  by 
the  positive  region  of  the  plant  nonlinearity.  D.F.G. 

1  B'  generates  the  factor  by  which  the  negative  step,  h.(T+) , 
must  be  multiplied.  The  comparator  selects  the  appropriate 
signal  during  the  correct  time  interval. 
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Fig.  7-2 

7-1-B  NEGATIVE  INPUTS 


For  negative  step  inputs  to  the  system,  the  initial 
input  to  the  plant  is  a  negative  step.  During  the  time 
interval,  T  <  t  <  2T,  a  positive  step  is  applied.  It 
can  be  seen  that  the  same  feedforward  nonlinear  elements 
(as  in  Fig.  7-2)  can  be  used  for  negative  inputs  as  long 
as  the  appropriate  scaling  factors  are  used.  The  scale 
factors  required  are  in  the  ratio  of  |h(0+)|  to  fh(T+)/  . 
Thus,  one  D.P.D.T.  relay  can  be  used  to  select  the 
appropriate  scaling  factors.  The  circuit  can  then  be 
modified  as  shown  in  Fig.  7-3. 


Fig.  7-3 
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7 -1-C  SUCCESSIVE  INPUTS 

If. successive  inputs  are  applied  to  the  system 
shown  in  Fig.  7-3,  deadbeat  response  will  not  be  obtained 
except  in  special  cases,  since  the  feedforward  network 
will  compensate  for  a  signal  which  is  equal  in  magnitude 
to  the  sum  of  all  the  previous  inputs. 

To  overcome  this  problem,  an  arrangement  must  be 
devised  by  which  the  feedforward  network  can  revert  to 
the  conditions  for  zero  input  after  the  system  has 
responded  to  the  previous  input.  If,  when  the  next  input 
is  applied,  the  feedforward  network  senses  only  this 
applied  input,  Then  the  system  will  respond  in  a  deadbeat 
fashion.  The  arrangement  used  to  do  this  uses  two  time 
delays  (see  Fig.  7-4).  The  system  as  shown  in  Fig.  7-4 
will  also  respond  with  minimal  response  when  any  of  the 
inputs  is  removed. 


Fig.  7-4 
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7-2  SYSTEM  WITH  'SANDWICHED'  NONLINEARITY 

With  reference  to  Fig.  4-5,  the  following  simplified 
system  will  be  considered  (Fig.  7-5). 


the  positive  portion  of  the  nonlinearity.  Thus  a  non- 
symmetrical  nonlinearity  affects  the  system  response 
only  when  negative  step  inputs  are  applied  or  when 
positive ’ inputs  are  removed.  To  compensate  for  these 
situations,  a  separate  set  of  nonlinear  functions  must 
be  derived.  Comparators  can  then  be  used  to  select  the 
appropriate  multiplying  factors  for  the  corresponding 
input  signals. 

By  using  two  time  delays  in  a  scheme  similar  to  the 
one  used  in  Fig.  7-4,  the  system  can  be  modified  to 
respond  to  successive  step  inputs  with  deadbeat  response. 
The  complete  system  is  shown  in  Fig.  7-6. 
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Fig .  7-6 

7-3  CONCLUSION 

In  this  chapter,  the  systems  were  compensated  so 
that  deadbeat  response  was  obtained  in  many  special 
situations.  The  compensation  used  here  may  or  may  not 
be  useful  depending  upon  the  inputs  being  considered 
and  the  response  required.  However,  using  the  methods 
outlined,  a  nonlinear  system  can  be  compensated  so 
that  the  response  is  minimal  for  almost  any  of  a 
variety  of  inputs.  The  only  restrictions  are  that: 

(1)  The  input  must  be  a  step  input. 

(2)  The  magnitude  of  the  step  must  be  less  than 
the  physical  limitations  of  the  system. 
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(3)  Successive  inputs  must  not  be  applied  within 
two  (or  in  some  cases  three)  sampling  periods  of  each 
other . 

Other  than  these  restrictions,  feedforward 
compensation  can  be  used  to  compensate  a  system  to 
which  either  positive  or  negative  successive  step 
inputs  are  applied. 
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APPENDIX  A 


SIMULATION  OF  THE  DIGITAL  CONTROLLER 

Through  most  of  this  thesis,  the  Digital  Controller 
was  simulated  using  a  Digital  Controller  Box  which  was 
recently  designed  and  built  by  the  department  especially 
for  this  purpose.  An  explanation  of  its  operation  using 
block  diagrams,  will  be  given  in  this  appendix. 

The  two  main  requirements  in  simulating  a  digital 
controller  are  (1)  a  device  for  delaying  a  signal  for  1, 
2,  ...  ,  n  sampling  periods,  if  D(z)  is  of  the  form: 

"1  -n 

aQ  +  ai  2  +.  .  .+  anz”n 

D(z)  =  — - — — 1  — - - - 

1  +  b-^z  ^  .  .  .+bnz~n  , 

and  (2)  a  device  to  hold  an  instantaneous  signal  through¬ 
out  one  sampling  period  (zero-order  hold), 

A  unit  time  delay  can  be  simulated  by  using  two. 
zero-order  holds.  If  the  second  zero-order  hold  samples 
and  holds  the  output  of  the  preceeding  one,  immediately 
before  the  preceeding  one  samples  and  holds,  then  the 
output  of  the  second  is  equal  to  the  output  of  the  first 
delayed  by  one  sampling  period.  This  is  valid  only  on 
the  assumption  that  the  input  is  described  by  its  values 
at  the  sampling  instants. 

By  following  a  unit  time  delay  with  an  additional 
zero-order  hold  which  samples  an  instant  before  the 
preceeding  one,  an  additional  time  delay  is  generated. 
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In  Fig.  A-l  a  block  diagram  schematically  outlines 
the  method  used  in  generating  these  time  delays. 


Fig.  A-l 
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Timing  pulses ,  generated  by  an  accurate  timing  clock, 
trigger  the  first  of  a  series  of  monostable  multivibrators 
as  illustrated  in  Fig.  A-l.  The  output  of  each  stage  is 
either  a  200  or  a  40  microsecond  pulse.  Each  stage  is 
always  triggered  by  the  lagging  edge  of  the  preceding 
pulse.  The  zero-order  holds  are  triggered  by  the  200 
microsecond  pulses,  the  last  zero-order  hold  being 
triggered  by  the  first  pulse  and  so  on.  The  zero-order 
hold  used  to  hold  the  output  of  the  digital  controller 
is  triggered  at  the  same  time  as  the  first  zero-order 
hold  in  the  train. 
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APPENDIX  B 


DIGITAL  COMPUTER  PROGRAM 
In  Chapter  4,  a  method  was  outlined  wherein  the 
digital  computer  was  used  to  calculate  the  points  on 
the  nonlinear  curves  used  in  the  feedforward  network. 

The  flow  sheet  for  the  computer  program  as  well  as  the 
program  itself  are  shown  in  this  appendix.  The  symbols 
used  are  defined  below. 

A(N)  -  ordinate  break  points  of  the  nonlinearity.  (See 
B(N)  -  abscissa  break  points  of  the  nonlinearity.  Fig.  4-4) 
T(I,J)  -  time  required  for  X2(t)  to  rise  to  B(J+1)  during 
the  test  where  X2(T)  =  B(I+1). 

TT(I,J)  -  Time  required  for  x^Ct)  to  die  down  from  B(I+1) 
to  B(I-J+1) 

S(I,J)  =  EXP(-T(I s J) ) 

SS(I,J)  =  EXP(-TT(I ,  J)  ) 

J 

H(I)  corresponds  to  the  value  h(0  )  required  for  X2(T) 
to  equal  B(I+1). 

HH(I)  corresponds  to  the  value  h(T+)  required  for  X£(t) 
to  die  down  from  B(I+1)  to  zero  during  the  time 
interval,  T  <  t  <  2T. 

ARA(I)  -  area  under  the  curve  of  x^(t)  between  0 <  t<^  T. 
ARB(I)  -  area  Under  the  curve  of  X2(t)  between  T<t<2T. 
SL(I)  -  slope  of  the  line  segment  of  the  plant  nonlinearity 
between  A(I),B(I)  and  A(I+1) ,B(I+1) . 


V 


a  xiijrcmA 


m-  'oc-.-'  i.  ht;.-,  ■  :■ 


■ 

o  ;•  inioq  :  ol  joIbj  o::  :.v:  .  a&w  tt'.Jijiqmro  for}  i  gib 

9dJ 

.  •  •  '  :  •  .  i’’  :  ■■  : 

••  '■  i  !  •  3  81 


398)  .  ^aixBsnilfron  s.<i3  1c  3:  riioq  alBSTd  s^Bnibto-^  (M)A 


(£-£  .  3  ^  .  ^bibssni  Irion  »ri3  Ho  aci.rioq  >U>io  seajC£f:d b  «*  (H)fl 

gnisub  (I+t)a  ob  98 z?  oJ  (3)$*  loi  bsxUi p»*  9ib±3  -  (LtI)T 

.  (X+‘  )d  **>  (T)CX  3" 9  *V7  3893  9fi3 


( X+I )a  moil  nvrob  s  b  o3  (3  ,x  *rbri  b^'iiu  .>97  smlT  -  (UfI)TT 


z 

* 


(I+L-I)a  03 

((l,I)T-)1Xa  »  (1,1)2 
((L<I)TT“)<IX3  *  (Lt 1)33 

;9  (I  )H 

•  * 

.(1+1)0  :p.upc.  03 

■ 

(3)^x  ioi  bs^lnpoi  (  ■') ii  swil  >v  sito  63  abnoqaesioo  (I)HH 
»c  orto  i  ril il  i  <nos  x3  (  1  /  mofi  nwob  3lrb  o3  \ 

’’  3  ri  .  i  :**»  -9  U •• : 


,f  ->3  >0  n99w39d  (3)*x  9visjo  e  rj  nobnu  B91B  -  (I) AHA 
S  *:•  T  iR>9>73ocf  f  ■  c  to  ovtuo  9d3  loj  ,rjtr  usxa  (I)fiHA 


^li-xsofri. Inon  3nBlq  orb  io  Ir^mgQa  ©nil  srto  iio  sqola  -  (I)J8 


.  (i+I)fit  (1+I)A  bna  (1)8,  (I) A  nssvrtsd 


-  62  - 


C ( I)  -  abscissa  intercept  of  the  line  segment  with  slope 
SL(I)  as  described  above. 

N  corresponds  to  the  number  of  line  segments  making  up 
the  plant  nonlinearity. 

AK  =  Kq 
BK  =  Kx 

RO  - 

o 

FAK  -  Rf  -  Feedforward  multiplying  constant  during 
the  first  sampling  period. 

FBK  =  Kf 
ER  -  e(T+) 

The  flow  diagram  for  the  computer  program  is  shown 
in  Fig.  B-l,  the  program  itself  being  shown  in  Fig.  B-2. 
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S JOB  345009 

$  IB JOB  FOR COM  "HODECK iNOGO 
$  I BFTC  GRK 

DIMENSION  A(15)*B(15)*T(15»15)(TT (15*15)* 

IS  C  15*15 ) *SS( 15*15 ) » H ( 1 4 ) * 

2HH  (  15)  .ARA(  14)  *  ARB  (  14)*SL(14)*C(14) 
READ(5*_10HA(  I  )♦  1  =  1*15) 

READ(5»10MB(  I  )  *  1  =  1*15) 

10  FORMAT  (6E12.5) 

N=  1 4 

AK=1. 58197671 
BK  =  1 »  3922 1 1 1 9 
UO_20„L=1  »N 

SL (  I  )  s ( B ( I +1 ) -B ( I ) ) / ( A ( 1  +  1 ) -A ( I )  ) 

C(  I  )  =  (A( 1+1 )*B( I )~B(  I  +  l )*A( I )  )/( A( I  +  l )-AC  I )  ) 
H(  I )=A( I  +  l ) / ( l.-EXF( -1 •  )  ) 

HH (  I  )  = A ( I + 1 ) ~H (  I ) 

WR I TE ( 6 ♦ 1 2  >  SL(I)*C<  I)*H(I)*HH(I  ) 

12..  .FORMAT  (  1X.4E16.6) 

ARA ( I ) =0. 

ARB (  I  ) =0. 

DO  30  J=1»I 

S(  I  »  J  )  =  (  H  (  I )-A( J  +  l )  ) / H  <  I  ) 

T ( I , J)=-ALOG(S( I »J) ) 

J I = I -J+l 

SS (  I » J )  =  ( H ( I ) -A (  I  +  l )+A< JI  )  ) / H ( I ) 

TT (  I » J ) =~ALOG( SS ( I » J )  ) 

WRITE(6*14)S(I*J)»T(I,J)*SS(I»J)»TT(I,J) 

14  FORMAT ( IX  *4E16.6 ) 

IF ( J.GT.l )  GO  TO  15 

ARA (  I  ) = AR A (  I  )  +H< I )* ( T ( I  * J ) +S  C  I » J ) -1 .  ) 

ARB(  I  )  = ARB (  I  )  +  (SL( I ) *HH ( 1 ) +C (  I )  ) * 

1 TT (  I  ,J)-SL(  I )*H( I )*(SS( 1  *  J ) -1 .  ) 

GO  TO  30 

15  AR A (  I  ) = AR A ( I )  +  ( SL ( J ) *H ( 1 ) +C ( J )  ) * 

1 ( T (  I » J ) -T {  I  * J-l )  )+SL ( J) *H(  I  ) * 

2 ( S (  I  ♦ J ) -S (  I , J-l )  ) 

ARB (  I  ) = ARB (  I  )  +  ( SL ( J I  ) *HH ( I ) +C ( J I )  ) * 

1 ( TT (  I  * J  )-TT(  I  * J-l )  ) -SL ( JI  ) *H (  I ) 

2*( SS(  I  * JJ-SSI I »J-1 )  ) 

30  CONTINUE 

WR I T E ( 6  >  3  5 )  ARA (  I ) * ARB (  I  ) 

35  FORMAT ( 1X,2E16.6 ) 

RO=ARA (I j+ARB (I ) 

FAK  =  H(  I  )  / ( RO*AK ) 

FBK=HH ( I ) / ( RO*BK ) 

ER=ARA( 1 ) /RO 

WRITE (6 *2  5  )  RO*FAK*FBK*ARA( I ) *ER 
25  FORMAT ( 1X*5E13*5) 

20  CONTINUE 
END 


Fig.  B-2 
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